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Numerous components operating in reusable space propulsion systems such as the 
Space Shuttle Main Engine (SSME) are exposed to high pressure gaseous hydrogen 
environments. Shown (Figure 1.0) is the SSME Powerhead. Flow areas and passages 
in the fuel turbopump, fuel and oxidizer preburners, main combustion chamber and 
injector assembly contain high pressure hydrogen either high in purity or as 
hydrogen rich steam. As can be seen, this includes many components such as turbine 
disks and blades, impellers, housings, ducts, etc. Accurate constitutive and 
damage material models applicable to high pressure hydrogen environments are 
therefore needed for engine design and analysis. Existing constitutive and cyclic 
crack initiation models have been evaluated only for conditions of oxidizing 
environments. The purpose of this effort (Figure 2.0) is to evaluate these models 
for applicability to high pressure hydrogen environments. 

Material behavior is known to be significantly affected by high pressure hydrogen 
environments (Figure 3.0). For the materials typically employed in rocket engine 
applications, hydrogen environment embrittlement is most severe near room 
temperature (Ref. 1). Tensile ductility and notched bar ultimate tensile strength 
reductions are typically reduced by the hydrogen environment. The room temperature 
embrittlement effects increase with increasing pressure (Ref. 2). Although the 
room temperature tensile properties of the superalloys are relatively insensitive 
to strain rate effects, hydrogen environment embrittlement is strain rate 
sensitive. Low cycle fatigue lives, especially crack initiation lives, are also 
reduced by the hydrogen environment (Ref. 3). Thermal activation and time 
dependent deformation (creep) also become important as temperature increases 
(Ref. 4). Indications are that the presence of hydrogen can accelerate creep 
effects by enhancing dislocation mobility (Ref. 5). 

The program flow chart (Figure 4.0) is shown for each material. Constitutive 
modeling effort precedes life modeling. Therefore constitutive test results can be 
utilized to provide a preliminary projection of crack initiation results, ductility 
normalized, (Ref. 6) as well as identify target test values for a, e, e, 

etc. During life testing, specimen constitutive response will be measured to 

provide an added database for constitutive model verification including any 
modifications. The experimental phase will begin with the selected isotropic 
material, followed by the anisotropic material. Both helium and hydrogen testing 
will be performed to obtain contrasting inert and aggressive environment test 
results . 

Alloy and material model selections as well as test conditions have been defined. 
Currently, existing data is being surveyed and compiled to provide initial 

estimates for model constants and to guide test details. Inconel 718 was selected 
as the isotropic material due to its extensive usage on the SSME, its 

susceptibility to Hydrogen Environment Embrittlement (HEE) and the existence of 
service-related hydrogen assisted cracking. PWA 1480 was selected as the single 
crystal material, due mainly to its development for potential usage on the SSME as 
a blade material and its known susceptibility to HEE. The nominal properties of 
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these two materials are presented (Figure 5.0) along with the test specimen 
geometry (Figure 6.0). The maximum test temperatures and pressur < Fl 9 ure 
correspond to maximum usage on the SSME. Due to the known enhanced HEE effects at 
room temperature ambient temperature testing shall also be performed. 

Recent emphasis on constitutive model development has concentrated on unified 
viscoplastic models wherein plastic, creep and relaxation strains are treated as 
one strain component, the inelastic strain. Constitutive model characterization 
tests (Figure 8.0) will be carried out at both room and elevated temperature in 
helium and hydrogen environments. The Walker model was chosen for the 1s ° tr °J2J 
material. Inconel 718 and Walker’s anisotropic model was selected _ for PWA 1480 
(Fiqure 9 0 Refs 7, 8, 9). A modification based on Valanis' material time scale 
has been 'developed in-house to extend these rate dependent models into the rate 
independent (low temperature) regime (Figure 10.0). 

Recent crack initiation model development has emphasized life prediction in 
oxidizing environments under such loading complexities as thermomechamcal fatigue, 
cyclic creep, multiaxial and variable amplitude loading. For assessment in a 
hydrogen environment, the total strain version of strain range partitioning 
(TS-SRP) (Refs. 10 f 11) was selected for constant amplitude loading conditions. 
Isothermal testing will be performed to determine the four basic inelastic strain- 
range life relations, PP, PC, CP, CC as well as the elastic strain-range life 
relation (Figure 11.0). Due to the preponderance of thermal gradient driven strain 
cvclinq of SSME hardware, the most important non-isothermal loop types to consider 
are out-of-phase PP and PC. To better simulate the impact of thermal mechanical 
fatique (TMF) , out-of-phase bithermal testing will be performed in PP cycling 
(Fiqure 12.0). Reasonable estimates of strain cycle history relative to TMF can be 
obtained where free thermal expansion strains can be easily subtracted from the 

total strain. 

Cumulative damage will be assessed via the advanced NA D S V e / V D e1 f ° Pe ?or <Je !h; 

Curve Approach (Figure 13.0) and Double Linear Oamage Rule (Ref. 12). The nature 
of these models is to accumulate damage nonlinearly with high-low load sequences 
being more damaging than low-high. Isothermal tests will be run to provide 
baseline data and a bithermal LCF test followed by HCF is planned as a verification 

experiment. 

The experimental program as proposed shall include uniaxial monotonic, cyclic (both 
isothermal and bithermal), creep and relaxation tests to invest! gate the 
applicability of these chosen models and to provide a basis for nece> a y 
modifications (Figure 14) arising from model deficiencies should they occur. 


*Work performed under NASA contract NAS3-26130 
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[ I CURT. 1.0 

COMPONENT S f XPOSf.l) 10 HIGH PRESSURE IIYDHOGLN 


• DETERMINE REGIMES OF MATERIAL BEHAVIOR ALTERED BY HIGH PRESSURE H 2 

• EVALUATE EXISTING MATERIAL MODELS FOR HIGH PRESSURE H 2 APPLICABILITY 

• PROPOSE MODEL MODIFICATIONS AS REQUIRED 




b/c 





CONSTITUTIVE MODEL LIFE MODEL CUMULATIVE DAMAGE MODEL 


FIGURE 2.0 

MATERIAL MODELING TASK OBJECTIVES 
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Ilydrogen/llelium Redoclion of Area 
for Inconel 718* 



Influence of Environmental 
Hydrogen on Inconel 718 Life 



TamfMratura ( C) 

*51 ?MPj HjOrojen 

HEE of Inconel 718 in Hydrogen 
as a Function uf Temperature* 



Creep Stress- Rupture of DS MaR- 
M246 in 760'C, 34.5 MPa Hydrogen and Helium 
Environments at 758.4 MPa Stress 


FIGURE 5.0 

MATERIAL ISSUES RflAiri) TO HIGH PRESSURE HYDROGEN 



FIGURE 4.0 PROGRAM FLOW (EACH MATERIAL) 
H2 AND HE TESTING 
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I SOI ROP [ C MATERIAL 
INCONEL 718 


SINGLE CRYSTAL MATERIAL 
PWA 1480 


NOMINAL COMPO S 1 T IQN_(^ LLGHJU£RC£MI1 


CR - 19 

MO 

- 3 

NB + T2 - 5 

TI 

- 1.0 

AL - 0.5 

FE 

- 18.5 

Nl - BAL. 




CR - 10 CO - 5 
W - 4 TI - 1.5 

TA - 12 AL - S 
N1 - BAL. 


HEAI_IREAJU£NI 


HEAT_JB£AItt£NJ _AMD_£RQ££SSI NG 

(TYPICAL TIMES AND TEMPERATURES) 


1900°f / 50 M1N/AIR COOL 
1400°F/L0 HRS/FURNACE COOL TO 
1200°r/HOU) TO 20 HOURS TOTAL 


2350°F/4 HOURS/AIR COOL 

HOT ISOSTATIC PRESS/4 HOURS/2350 0 F /I 5 KSI 

1975°F/4 HOURS 

1600°F/32 HOURS 


TENSILE PROPERTIES (TYPICAL), KSI 

RL 65Q •£ 

HE H 2 HE 

F jy 160 160 140 

F tu 185 185 160 

ELONG . X 20 11 20 


TENSILE PROPERTIES (TYPICAL) <001>, KSI 
RT m°JL 

HE H2 HE 
F ty 145 140 100 

F tu 165 150 140 

ELONG. % 5.4 2.5 12 


FIGURE 5.0 
MATERIAL SELECT IONS 


SPECIMEN DIMENSIONS 


HTFV Bultcm Head Fatigue Sample 


3.50 



7 

0 . 675 <#> 
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|| VESSEL DESCR1PI ION/ 1 NS I ROME NT AT ION DESCRIPTION 

• 10.000 P$I/1800°F FATIGUE SYSTEM 

• INTERNAL LOAD AND ON-SPECIMEN STRAIN CAPABILITY 

• SELF COMPENSATING (TEMP.. PRESSURE, ENVIRONMENT) CAPACITANCE EXTENS0ME1LR 

• FULL COMPUTER TEST CONTROL AND DATA MANAGEMENT 


FIGURE 6.0 

TEST SPECIMEN AND VESSEL DESCRIPTION 
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ISOTROPIC MAT! RIAL - INCONEL 718 


SINGLE CRYSTAL MATERIAL 


• MAXIMUM USAGE (SSMEI 

• TURBINE HOUSING AND OIJ1IU MANIFOLD 

• PRESSURE = 2N-3'l MPa 

• HOI GAS TEMP = fsSO 0 ] C 

• PUMP INLET TO OIJ1LET 

• PRESSURE = [SO”! MPA 

• IEMP = 20° C 

• MINIMUM TEMP 

• l2jr] C 


PWA IN 80 

• MAXIMUM USAGE (SSME) 

• HIGH PRESSURE TURBOPUMP BLADES 

• PRESSURE =[3N;>1 MPA 

• TEMP =[870*] C 


• MINIMUM TEMP 
• [ 20®] C 




I IGIJRL 7.0 

If SI IIMPI RAIURI S AN!) PRfSSHRIS 


Monotonic Tensile Tests 



Cyclic Tests 




Control History 



A 

V 


Creep/Relaxation Tests 



•i <» 


o' 


Figure 8.0 

CONSTITUTIVE MODEL CHARACTERIZATION TESTS 
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ISOTRQ P! CLCQNSTIMiTI YE JV1Q DEL 


S I NG LE CR YSTAL CON STITUTI V E MODEL 


jU a tonal Selector!: INCO-718 
Mod el S elected: Walker’s Isotropic Model 
Mollification: Employ Material Time Concept 
Formulation: 

* d . tlz 

n= I: u-. «"<' i = s - 


Material S elected: PWA-148U 

Model Selected: Walker’s Single Crystal Model 
Modification: Employ Material Time Concept 


Formula tio n: 


• The Flow Rule; 


) w-*i 


• The Evolution Equations: 

K = K ] - 

* . . ,r - fln, < 1 

a, - (».. tn,) •„ 



0=0 


dz 


S^°j) 


■ The Hooke’s Law: 

<7.j = («.; - «fj) ■» 

» -V (.„ - I' i) t- ' V 


• The FlowRule; on llie r‘* slip system of the 1 2 oct.ilinlra! and 
(j «■ ubit. systems ( w r below is the Schmid shear stress on the slip plane) 

of) 

• The Evolution Equations for u r and_A' r 

w r = p X 7r ~Pl |7 f |<*\ - i <->r P Wr 

Ar - |£l;% t (1 -q)6,*) -v(Kr- A‘)||7,|J - M hf - h-y 

• The Inelastic Strain Rate: 

t3 « 

«' = EX I* + XX 7, 

• The Hooke’s Law: 

R.,u|#u - 


Figure 9.0 WALKER’S UNIFIED CONSTITUTIVE MODELS 


Basicjdeal 

• Material experiences deformation and remembers its deformation history. 

• Tims, the time felt by material should also include the effect of deformation. 

Expressio n of the M ateria l Time Increment 

dz = y/Kjkt<-ij<kl + <jdt 2 
where A,^ and <y are material parameters. 

Proposed Use of the Mat erial Tim e 


Change the time derivatives in a viscoplastio formulation into derivatives with respect 
to the material time. 

* It reduces lo the original viscoplastic formulation when A tJ *| is null and <j = 1 (used for high 
temperature cases) 

4 It has no rale dependence while plasticity remains when \ tJ ki is constant and ry = 0 (applicable 
for low temperatures) 


Figure 10.0 THE MATERIAL TIME CONCEPT 
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INELASTIC SRP CYCLE 
TYPES/IDEALIZED LOOPS 
(Re = -1.0) 



lit PP TYP{ CYCtC. 

smtssi 



CSttP 

fc» pc rvpt cycle 


stwes s 

f CREEP 



(bicp type cycle 


smESs 
I CRtf P 



CREEP 

Idt CC rrPf CYCLE. 


TYPICAL PARTITIONED STRAIN 
RANGE ~ L 1 EE RELATIONS 



INTERACTION DAMAGE RULE: 

F. . , 

l jr 1 = 1 . ij = PP, pc, cp. cc 

i j to 


FIGURE 11.0 

DAMAGE MODEL TS-SRP ( ISOTHERMAL) 
Rc = -1.0 AND Rc = 0.0 (PP) 


OUT-OF-PHASE TMF 


OUT-OF-PHASE BITHERMAL 


TEMPERATURE 
AND STRAIN 
WAVEFORMS 




HYSTERESIS 

LOOPS 


STRESS 



STRAIN 



FIGURE 12.0 

COMPARISON BETWEEN HI T Ilf RMAI AND IMF CYCI INC, 
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TWO STEP LOADING PATTERN 



FIGURE; 13.0 
DAMAGF. CURVE: APPROACH 
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